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Abstract 

An analytical model that represents the collapse of a massless 
scalar wave packet with continuous self-similarity is constructed, and 
critical phenomena are found. In the supercritical case, the mass of 
black holes is finite and has the form M oc (p — p*) 7 , with 7 = 1/2. 
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I. Introduction 

Recently, Choptuik |l] studied the collapse of a massless scalar field with 
spherical symmetry, and found the following intriguing features: Let the 
initial distribution of the massless scalar field be parametrized smoothly by 
a parameter p that characterizes the strength of the initial conditions, such 
that the collapse of the scalar field with the initial data p > p* forms a 
black hole, while the one with p < p* does not. Then, it was found that: 
(i) the critical solution with p = p* is universal in the sense that in all the 
one-parameter families of the solutions considered it approaches an identical 
spacetime; (ii) the critical solution has a discrete self-similarity (DSS); (iii) 
near the critical solution (but with p > p*), the black hole mass is given by 

M BH = K(p-p*)\ 

where K is a family-dependent constant, but 7 is an apparently universal 
scaling exponent, which has been numerically determined as 7 ~ 0.37. These 
phenomena were soon found also in the collapse of axisymmetric gravitational 
waves [0 as well as in the one of radiation fluid ||. Therefore, it seems that 
the phenomena are not due to the particular choice of the matter fields but 
rather are generic features of General Relativity. Further numerical evidences 
to support this conclusion are given in Q. 

Parallelly to the above numerical investigations, there have been analyt- 

2 



ical efforts to understand the physics behind these phenomena || ||, [7J ^ 
H |lOfl . While the universality of the critical solution and its self-similarity 
(echoing) have been found in most cases considered, there exists dispute 
over the universal scaling exponent 7. Maison || showed that 7 is matter- 
dependent. For the collapse of the perfect fluid with the equation of state 
P = kp, it strongly depends on k, where P and p are respectively the pres- 
sure and energy density of the fluid, and A; is a constant. This dependence is 
further shown in JK| . The same conclusion was also obtained numerically by 



Eardley and co-workers [IT|. Thus, one might expect that 7 is universal only 
within a particular family of matter fields. However, even in this sense the 
analytical results are not consistent with the numerical ones. In particular, 
it was shown |5|, |6| that for the massless scalar field the exponent 7 is 0.5, 
instead of its numerical value 7 ~ 0.37. Moreover, Maison || showed that 
the value k = 0.88 is the maximal one for which a regular self-similar solu- 
tion could be found for the perfect fluid. As we know, a massless scalar field 
is energetically equivalent to a perfect fluid with k = 1 |12| , thus, Maison's 
results seem also in conflict with the ones obtained in || |6fl. 

It should be noted that the model considered in || ||, which will be re- 
ferred to as the BONT model, is not asymptotically flat, and the exponent 7 
is obtained in a rather unusual way. Moreover, except for the one of Gund- 
lach H, all the analytic models, including the BONT model, have continuous 
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self-similarity (CSS), in contrast to the original model of Choptuik that has 
DSS 0. Thus, it is not clear whether the difference in 7 obtained above is 
due to the different self-similarities or due to the non-asymptotic flatness of 
the spacetime. If it is because of the former, we will have a new classifica- 
tion to the critical behavior, and many new solutions that represent critical 
phenomena with different similarities need to be found. To resolve the above 
problem, in this paper we shall make a "surgery" to the BONT model [0]. 
That is, we shall cut the spacetime along a null hypersurface and then join it 
to an out-going Vaidya solution, so the resulted spacetime is asymptotically 
flat and has a finite mass. Clearly, such obtained solution will represent the 
collapse of a massless scalar wave packet, which is radiating as it collapses. 
This model is more realistic and more comparable with the one of Choptuik 
but with different self- similarities. Specifically, the paper is organized as fol- 
lows: In Sec. II, we briefly review the main properties of the BONT model, 
while in Sec. Ill we cut the BONT spacetime along a null hypersurface and 
then join it with an out-going Vaidya solution. In Sec. IV we consider a 
particular case of the solutions obtained in Sec. Ill, which represents the 
collapse of a scalar wave packet. Critical phenomena are found with the 
exponent 7 = 1/2. Finally, the paper is closed by Sec. V, where our main 
conclusions are presented. 
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II. The BONT model 



The BONT model is described by the solutions [] 



ds 2 = —G(u, v)dudv + r 2 (u, v)d 2 Q, 



(1) 



where d 2 Q = d9 2 + sin 



2 QA,r2 



and the metric coefficients are given by 



r(u, v) 
G{u,v) 



u 2 - 2uv + Ab 2 v 2 ] V2 H(v) + ^ [2{a(v) - a(0)) - u] #(-«), 



#(v) + 2a'(v)iZ'(-v), 



(2) 



where b 2 is an arbitrary constant, and H(x) denotes the Heaviside function, 
which is one for x > and zero for x < 0. A prime denotes the ordinary 
differentiation, and a(v) is an arbitrary function subject to a'{v) > and 
a'(0) = 1/2. One can show that the hypersurface v — is free of any matter 
and represents a boundary surface ||Plfl . The corresponding massless scalar 
field is given by 



1 , 

±— = In 
V2 



(u — v ) — s/l — Ab 2 v 



H(v). 



(3) 



(u-v) + VI - 4b 2 v 
Note the slight difference in the notations used here and the ones used in || §]. 
From the above expressions one can show that the spacetime is Minkowski 

in the region u < 0,v < 0, while in the region w < 0, f > it represents a 

1 This class of solutions was first found in J13| , but their physical interpretation of 
representing critical phenomena was first given in |5|, ^) . The Roberts solutions are contin- 
uous self-similar, since they possess the Killing vector, X — udu + vdv with the property 
Lx9fj.ii = ^9^u, where Lx denotes the Lie derivative with respect to X. 
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collapsing massless scalar wave. When 62 < 0, the scalar wave collapses into 
a spacetime singularity on the hypersurface u = —[VI — 46 2 — l]i>, which 
is preceded by an apparent horizon at u = ^o 2 v - Thus, the corresponding 
solutions represent the formation of black holes and are supercritical. When 
62 = 0, the singularity coincides with the apparent horizon on u = and 
becomes null. This solution is critical, which separates the supercritical 
solutions from the subcritical ones. The subcritical solutions are these with 
< 62 < 1/4. It can be shown that in the latter case the scalar wave first 
collapses and then disperses into infinity, without forming black holes but 
leaving behind a Minkowski region u, v > 0, in which the metric takes the 
form of Eq.([I]) with 



where b(u) is an arbitrary function, subject to b'{u) > 0,6'(0) = 1/(4-^/62). 

One can show that the hypersurface u = 0, v > is also a boundary surface. 
In the region u < 0, v > 0, the local mass of the scalar field is given by 



where r is given by Eq.(|2|). Clearly, on the apparent horizon u = 4&2^ the 
mass becomes unbounded as v — > +00 for 6 2 < 0. That is, the spacetime 
fails to be asymptotically flat for the supercritical case. As a result, the total 









6 



mass of the black hole can not be written in a power-law form in terms of 
initial data, whereby the exponent 7 can be read out. For more details, we 
refer the readers to @, |§] . 

III. Matching the Roberts solutions to the out- 
going Vaidya dust solution 

To circle the above problem, we shall cut the spacetime along the hypersur- 
face, say, v = v > 0, and keep the region v < v , while the region v > v will 



be replaced by that of out-going dust Vaidya solution |15| . To do the match- 
ing, following Barrabes and Israel [fRJ, we first write the Roberts solutions 
in the form 

ds 2 _ = -e^-dv{f.e^-dv - 2dr) + r 2 d 2 Q, (6) 

where v is the Eddington advanced time, and for Roberts' solutions @ in 
the region u < 0, v > 0, the metric coefficients are given by 

f _ = [(1 - 4 t2 y + 1 [(1 _ ihy + 4r2] m _ (1 _ 46g)p | _ 

e>- = — (7) 

[(1 - Ab 2 )v 2 + 4r 2 ] 1/2 

Now we restrict Eq.(0) valid only in the region u < 0,0 < v < vq. Then, 
the normal to the hypersurface $~ = v — Vq = is n~ = a^d^^ = a -1 ^, 
where a is a negative function. From n~ we can introduce a "transverse" 
null vector iV~ by requiring N^N~ X = 0, and N^n~ x = —1. Without loss 



of generality, we assume that N~ takes the form N~ = N~5^ + N~5 r ^, and 
choose the arbitrary function a as a = —exp{—ip-}. Then, it is easy to show 
that N~ is given by 

n- = -^'-6; + 6;. (8) 

Choosing the coordinates r, 9, and if as the three intrinsic coordinates 
£ a = (r, 9, (p), (a = 1, 2, 3) on the hypersurface v = t>o, we find 

e (l) = ^f) e (2) = e (3) = (9) 

where eT^ = dx^L/d!; a . Then, it can be shown that the "transverse" extrinsic 
curvature, defined by [[16] 

K ab = -N,e\ b) (v„ef a) ) , (10) 

takes the form 

{dip T f T f I 
— dr^' ~2~' ~ sin29 \ ' ^ 

Note that in calculating the above equation, we have not used the particular 
expressions ([7p for the functions ^_ and /_. Thus, it is valid for the general 
case. 

On the other hand, the out-going Vaidya solution ]T5| can be written in 
the form 

ds\ = -e^+dU{f + e^+dU + 2dr) + r 2 d 2 Q, (12) 



where 

/ + = 1-5=22, , )+ = , (13) 
r 

and U is the Eddington retarded time, which is in general the function of 
u appearing in Eq.(l), and m(U) is the local mass of the out-going Vaidya 
dust. The corresponding energy-momentum tensor is given by 



r 2 dU K ' 



In the following, we shall take metric ( O ) as valid only in the region u < 
0, u > Vq. To have our results more applicable, for the moment we shall not 
restrict ourselves to the particular solution ([TBI). The hypersurface v = v in 
the (Z7, r)-coordinates can be written as $ + = U — Uq(t) = 0, where Uq(t) is 
a solution of the equation 



dU 2 



— e"^, (« = «„). (15) 



Then, the normal to the surface is given by 



where (3 is a negative otherwise arbitrary function. From n+ we can also 
introduce the "transverse" null vector ]V+ by requiring N^N +X = 0, and 
N^n +X = — 1. It can be shown that it takes the form 

N+ = ^e 2 ^<. (16) 



On the other hand, we also have 



e Oj; -t- o r , e/ 2 ) — o e , e 



(17) 



°(1) — u C/^"r> °(2) — U 0J °(3) — *V 

where et^ = dx+/d£ a . To be sure that the two "transverse" vectors 
defined in the two faces of the hypersurface v = v$ represent the same vector, 
we need to impose the condition 



iV A e (a) 



V=VQ 



which requires that the function /3 has to be (3 = — exp{— ip + }. Once and 
et* are given, using Eq. (|It]) we can calculate the corresponding "transverse" 
extrinsic curvature, which in the present case takes the form 



K 



ab ^ dia 9- 



2e-*+df + rf + rf. 



fl dU' 2 ' 2 
Then, from Eqs. (|TT|) and (Hp, we find that 



±sin 2 fl . 



(18) 



lab 



2 (Kb - Kb) 
2dtp_ 4e"^+ 9/ 4 



<9r 



W + K/+ - /-)(« + sm 2 0<W). (19) 



Once 7 a b is given, using the formula [TIJ 

1 



r 



06 



-s 



ab 



167T 



(<£ c ^ d + g?l a l c - g?l c l d - g?l a l b )icd, 



(20) 



we can calculate the surface energy-momentum tensor r ab on the null hyper- 
surface v — Vq, which now can be written as 



T ab = al a l b + Pg 



a lb 



ab 
/* 5 



(21) 
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where 



and 



a 



/+-/- 



57rr 



P = — 



1 (2e-^+dj+ d%l). 



in 



fl dU dr 



(22) 



(SX + sin" 2 (9^) , 



r 



i% = o. 



(23) 



The function a represents the surface energy density of the null shell, and 



P the pressures in the 9- and (^-directions. Note that in [16| the case of a 
null shell was also considered. But, there they used the same null coordinate 
in both sides of the shell, while here we use the retarded null coordinate in 
one side of the shell and the advanced null coordinate in the other side [cf. 
Eqs.® and (0)]. 

For the particular solutions given by Eqs.(^) and (|13|), Eq. (|22"l) yields 



a 



Anr 2 



M(r) 



(1 - Ab 2 )v 



1 _ 4b 2 )vi + 4r 2 



1 idM(r) (1 - Ab 2 )vl{r - 2M) 

4yr(r - 2M) \ lTr 2r[(l - Ab 2 )vl + 4r 2 ] 



where 



M(r) = m(U)\ u=Uo(r) 
and Uq{t) is a solution of Eq. (p~5|) . 



(24) 



(25) 
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IV. Critical phenomena for the case P = 



To study the general shell given by Eq. (|24j) , it is found very complicate. In 



this section, we shall consider the case where P = 0, i.e., 



dM(r) (1 - 4b 2 )v%(r - 2M) 



(26) 



dr 2r[(l -46 2 )fo +4r 2 ]' 
Integrating the above equation, we find that 

M(r) = - \p\jAp* 2 + r 2 - 2p* 2 J , (27) 

where p is the integration constant, and 



V = 1 v . (28) 

At the past null infinity (r — > +oo), Eq.(^) shows that 

M(r -> +oo) = p. (29) 

That is, the parameter p represents the total initial mass of the massless 
scalar wave packet and the null shell, with which they collapse. As r — > + , 
M(r) has the following asymptotic behavior 

r +oo, p>p*, 

M(r) -> <^ 0, p = p*, (30) 
[ — oo, p < p*. 

On the other hand, it is well-known that the apparent horizon at r = 2M(r) 
of the out-going Vaidya solution always coincides with its future event hori- 
zon. Thus, by comparing the mass M(r) with r/2 we can tell whether the 
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collapse forms a black hole or not, 

M(r)-- = \ - 9r AP + r [ 4 (p 2 -j9* 2 )-r 2 l. (31) 

v; 2 \4r 2 [p+ (V 2 +r 2 )] J L VF F y J v ; 

Clearly, only when p > p*, the scalar field and the null shell will collapse 
inside the event horizon at 



r AH = 2y/p> -p* 2 . (32) 

When p = p*, M(r) = r/2 is possible only at the origin, r = 0, where a 
zero-mass singularity is formed. Thus, the solution with p = p* represents 
the critical solution that separates the supercritical solutions (p > p*) from 
the subcritical ones (p < p*). In the subcritical case, M{r) is always greater 
than r/2, and the collapse never forms a black hole [cf. Fig.l]. In the 
latter case, the region u, v > should be replaced by the Minkowski solution 
(U). As shown in the last section, the matching across the hypersurface 
u — 0, < v < Vq is smoothy, i.e., no matter appears on it. To show that 
it is also the case on the hypersurface u = 0, v > Vq, which separates the 
Vaidya solution (|T3| ) from the Minkowski (|4]), we first make the coordinate 
transformation U = U (it), and then write the metric (O) in terms of u. Using 



the results obtained in [16], one can show that to have a smooth matching 



we have to impose the condition 



U'(0) = -±=, M(r) | u= o = 
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Clearly, by properly choosing the function-dependence of U on u, the first 
condition can be always satisfied. On the other hand, from Eqs. (||]) and 
one can show that the last condition is also satisfied identically. Therefore, 
the matching of the Vaidya solution to the Minkowski one across the hy- 
persurface u — 0, v > v o is always possible for p < p*. The corresponding 
Penrose diagram for each of the three cases are shown in Fig. 2. 

On the apparent horizon r = tah, the total mass of the scalar wave packet 
and the shell is 

M BH = r -f = K(p-p*) 1/2 , (33) 
where K = (p + p*) 1 ^ 2 . The above expression shows that the black hole 
mass takes a power-law form with its exponent 7 being equal to 0.5. This 
is exactly the same value as that obtained in |5], ||. This result seems a 
little bit surprising, since our model is different from the one of [[5], In 
particular, in our case the collapse consists of two parts, one is the massless 
scalar wave packet and the other is the null shell, while in 0, |6| only a 
collapsing scalar wave exists. However, a more detailed investigation shows 
that the existence of the shell does not affect our main conclusions. In fact, 
by properly choosing the parameters in our model the shell can be made 
completely disappear, without affecting the main properties of the collapse. 
To show this, let us write the total mass M(r) as 

M(r) = M (j> + M shelu (34) 
14 



where denotes the total mass of the scalar wave packet, and M s h e ii the 
total mass of the shell, given respectively by 



M = ^l(2yV 2 +r*-v } 



v r 

v p- V 2 (Ap* 2 +r 2 \ 1/2 

M she u = < 5 > . (35) 

v { r 2 J 

Clearly, if we choose the parameters such that 

4p* 2 

v = , (36) 

V 

the null shell disappears. Moreover, using Eqs. (|28|) and (p6|), we find 

Ap* 



p — p 



-b 2 ). (37) 



p + p* 

Therefore, when the null shell is absent, the massless scalar wave packet will 
collapse to form a black hole for b 2 < 0, and first collapses and then disperses 
to infinity, without forming a black hole for 1/4 > b 2 > 0, [| . 

V. Concluding remarks 

In this paper, by "cut-paste" method, we have constructed a physically more 
realistic model for the collapse of a massless scalar wave packet, which is 
usually accompanied by a collapsing null dust shell. Critical phenomena 
are also found in these solutions, and the mass of the black holes is finite 
and takes the form M oc (p — p*) 1 ' . In our model, the exponent 7 is 0.5, 
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which is the same as that found in |5], but different from the numerical 
one of Choptuik, which is 0.37. However, the results do not contradict and 
rather show the fact that 7 is not only matter-dependent but also symmetry- 
dependent. In Choptuik's model, the critical solution has DSS, while in || || 
and ours it has CSS. 

On the other hand, quite recently Choptuik, Chmaj, and Bizoh |J7] stud- 
ied the collapse of a Yang-Mills field, and found two distinct critical solutions 
at the threshold of black hole formation. In one case, the critical solution 
has DSS, and the mass of black holes for the supercritical solutions takes 
a power-law form with the exponent 7 ~ 0.20. This class of solutions was 
referred to as Type II solutions. In the other case, the critical solution is the 
static n = 1 Bartnik-Mckinnon sphaleron [fl8 1 , and the formation of black 
holes always turns on at finite mass. This class of solutions was referred to 
as Type I solutions. 

Summarizing all the above results, the following seems to emerge: The 
solutions of gravitational collapse can be at least divided into three different 
types, a) Type II. A solutions. In these solutions, the mass of black holes 
always takes the form M oc (p — p*) lA . b) Type II. B solutions. In this class of 
solutions, the black hole mass also takes the power-law form, M oc (p — p*) lB , 
but the only difference to Type II. A is that the critical solution of Type II. A 
has DSS, while the one of Type II. B has CSS. Because of this difference, the 
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exponent 7 is in general also different, c) Type I solutions. These solutions 
have no self-similarities, neither DSS nor CSS, and the formation of black 
holes always turns on at finite mass. If the above classification is universal, 
a natural question is that: Do these three types of solutions saturate all the 
solutions that represent the formation of black holes in gravitational collapse? 
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Figure Captions 

Fig. 1 The qualitative behavior of the total mass M(r) defined by Eq.(|27j) 
in the text for the three different cases: (a) The supercritical case (p > p*); 
(b) the critical case (p = p*); and (c) the subcritical case (p < p*). 

Fig. 2 The Penrose diagram of the solutions with P = for the three 
different cases: (a) The supercritical case; (b) the critical case; and (c) the 
subcritical case. Double lines represent spacetime singularities. 
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